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Dynamics of a light quark in the field of static source (heavy-light meson) is studied using the 
nonlinear Dirac equation, derived recently. Special attention is paid to the contribution of magnetic 
correlators and it is found that it yields a significant increase of string tension at intermediate 
distances. The spectrum of heavy-light mesons is computed with account of this contribution and 
compared to experimental and lattice data. 
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I. INTRODUCTION 



The nature of the QCD string between static charges (the static string) was studied extensively both analytically 
OA ' JjHl ano - on the lattice ||-^|. It was shown in these papers that the static string is predominantly electric (the 
connected probing plaquette is used for the analysis) and the electric field is directed along the string axis. 

In terms of the Field Correlator Formalism (FCM) J7S the static string is made of the correlators of electric fields 
only, and recent analysis in terms of Casimir scaling Hltl shows that up to 1% accuracy the string is formed by the 
, bilocal correlator of electric fields D¥ (x, y) = (jj-trE\\ (x)&(x, y)E\\ (y)) (where $>(x, y) is the parallel transporter and 
E\\{x) is the component of electric field parallel to the string axis). Thus the confinement dynamics for heavy (static) 
quarks is defined by D®(x,y) with this accuracy. 

This fact was used to construct the effective Lagrangian fo r lig ht quarks in the field of the static charge (the heavy- 
light quark-antiquark system) . The analysis made in fllj] has shown that using this Lagrangian one can derive 
in the large N c limit the nonlinear and nonlocal Dirac equation for the light quark Green's function 

^ (-id-im-iM)S = I (1) 

D , where the kernel M = M(x,y) is proportional to the bilocal (Gaussian) field correlator 



Q 2 

(j^trF^(x)<P(x,y)F Xa (y)) = D^ M {x,y) 



and to the Green's function S. 

It was shown in |ll]|l^] that the scalar confinement occurs at large distances in the selfconsistent solution of (Q), 
signalling Chiral Symmetry Breaking(CSB). The subsequent analysis in |l5( ] has supported this conclusion and the 
spectrum of the heavy-light meson states was calculated together with first estimates of the chiral condensate. 

In |ll| it was assumed that the electric correlators are dominant in the string and the magnetic part could be 
neglected. On the other hand the analysis of the heavy quark mass case (i.e. of equation (|l]) with the replacement 
of M( D, S) b y M(D,So) where So is the free Green's function for the heavy mass m) done in jy], Appendix 5, 
and in 0,[15|, has shown that the magnetic correlators can also significantly contribute (at least in the regime when 
mTg <C 1, where T g is the slope of -D^act)- It is a purpose of the present paper to study systematically the role of 
magnetic correlators for the light quark mass case, m <C ^/a, and to make quantitative analysis in this case. As a 
byproduct of our analysis the case of heavy quark mass is reconsidered and some refinements of previous results are 
obtained. 



II. MAGNETIC FIELD CONTRIBUTION TO THE CONFINING STRING 

We study in this section the solution of Eq. ([[]) with the help of relativistic WKB approach, similarly as in [^l) . 
The kernel M (where only the bilocal field correlator is kept and the Gaussian form for it is assumed) can be written 
in the form 
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iM(h,x,y) =f IJ ,S(x,y)f 1/ J IJ , 1/ (x,y) 



(2) 



where 



and 



J^(x,y) = exp(-h 2 /4Tg)J^ (x,y), /i = x 4 - y 4 , (3) 



J 44 = J (B) (x,y) = xy/ E (x,y)^, (4) 



^fc = 4° = (xy<5 ife - y l x k )f M (x, y)-^=^- (5) 



Finally, 



/ B( M)(x,y) = [ ds [ dt{st) a 
Jo Jo 



(xs - yt) 



ex P( 4j| ) 



where a = for and 1 for Jm- 

In what follows only asymptotic values of Je,m will be of importance, with |x — y| <C |x|, |y| 3> T g , in which case 
one has 

/ B (x,x) = 3/M(x,x)-^p. (7) 

It should be noted at this point, that sub - and superscripts E and M refer to the correlators of color-electric and 
color-magnetic fields respectively. Due to the structure of Eqs. (||-6) for x = y — > 00 the kernel is dominated by the 
color-electric field contribution. On the basis of this, in (k| the magnetic part of M was disregarded. However, in 
what follows we will show that for light quarks the magnetic part plays an important role at intermediate distances 
and as a consequence it can affect strongly the lower lying states of the heavy-light systems. 

The kernel M, Eq. (||) , contains the light quark Green's function S 1 , which is a selfconsistent solution of Eq. (Q) . 
Following Rcf. pd| we can determine S assuming that M in the first approximation is instantaneous. In this case S 
has the spectral decomposition in terms of eigenvalues e n and cigenfunctions i/> n , 

S(h,x,y) = i|^e- £ «^(/ l )V„(x)V+(y) - 

I n+ 

£ e-^ h 6(~h)Mx)^(y) ) 74 = i{0(h)S™ - 9(-h)S^}^ (8) 



e n + = £n and e n - = —el where J2n+ ano - J2n- re f er to sums over positive and negative eigenvalues respectively. 
Insertion of (8) in (||) then yields for M 

M(/i,x,y) = -ij^SjvJ^ = 

9(h)[S^ l4 J^ - jff^^+W] - d{-h)[S(-h±J (E) - 4" ) 7iS ( - ) 747fc]. (9) 

To find the properties of M we replace S by Su n in Eq. (Q), obtaining in this way Mu n , where Su n is the quark 
Green's function of the linear Dirac equation ([!]) with the kernel M replaced by clx^ 4 ) (x — y) (i.e. the static Dirac 
equation with linear potential ax). We now demonstrate that Mn n indeed tends to <r|x| at large distances, and thus 
yields the a posteriori proof that the large distance dynamics of the heavy light quark system is governed by the linear 
static local confining potential. At the same time, in the framework of the same formalism, it will be shown that an 
intermediate distance a region appears, where the dynamics is again local in time and static but with a larger string 
tension due to the contribution of magnetic terms. For the spherical spinor expansion has the form 

s&(h,x,y) = ]T e - £ -Vn(x)VvUy) = 

(n±) 
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V e _L ( G n {x)G* n {y)Sl jlm tl* lm , -iG n (x)F*(y)n jlm n* Vm \ 

™, ;/-..i.r:C- :,/:,() ,..<>', . F„(.^F*(V)a,-,,„.n*,,... 1 UJ 



\iF n (x)G* n (y)n 3 i, m tt* lml F n (x)F*{y)£l jVm Q.* jVl 

Similarly as in Ref. JO], we may carry out the summation over partial waves in (|l0|) using the WKB approximation 
for the solutions G n ,F n . As exploited in pllpTJ the results for can simply be obtained from using the 

symmetry of e n + and e„- solutions, namely e\i— eii ^ = e n and (e„, G„, F n , re) <-> (— e n , F n ,G n , —re). We quote the 
final result of the WKB analysis for when |x — y| <§; |x| 

sW = ^(1-cob*)/> ±A °' J A ) (11) 

where A = (to + <rx)|/i|. The matrix elements X,X contribute a nongrowing part to M and will be of no interest to 
us in what follows, while A 1; A are defined as 



, 'dre-^ C -^^^ = ^K (V^^)e^Ma) (12) 

7T Jl Vr 2 - 1 7T 



and 



cos(aVr 2 - 1) = -e A = A i(«)- (13) 

^ Ji v t - 1 77 V A 2 + a 2 

Here a = (era; + m)\x — y| for ex 2 ^> 1 and |x — y| -C x. From expressions ([l2|) and (|l^) we see that Ao, Ai are 
normalized as 

/>oc />OC 

/ A Q {a)da= / Ai(a)do=l (14) 
Jo Jo 

and hence diagonal elements of behave as smeared 6- functions: 

|5W(fc,x,y)rf , y = ie- A ( 1±1 1t1 ). (15) 
Indeed for large era; the functions Ao, Ai decrease exponentially fast when \x — y\ increases. Consider now Eq. (Q), 

(-i'r M d M -im)S(h,x,y)-i [ e'^T dh'e-^ x+m ^ h - h '\x 



$ w (a:-tf). (16) 



In Eq. ( |16[ ) the integration J Af(x, z)S*(z, y)d 3 z has been carried out using Eq. (|15[) - The obtained equation is a 
Dirac equation with a time-dependent interaction, localized in configurational space. We may now take into account 
that at large x ^S> T g we have 

^ ) (x,x) = ^ r ; ^ = i™ (17) 

There exist two regions in x, where Eq. ( |l6| ) can be simplified further. Considering the integration over d/i' in (|l|), 
in the situation when era; 2 3> 1 and T g — > 0, we have the two possibilities 

i) (to + cra;)r 9 < 1 (18) 
m) (m + CTx)T s >l. (19) 
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In the first case, i.e. when (|T^) holds, the leading contribution to the integral over dh' come from the region \h' — h\ < 
2T g because of the first factor in the integrand of Eq. (|l6|). Since the remaining factors vary smoothly over this 
region, provided (|||) is satisfied, we may replace h' = h in these factors, including S(h',x.,y). In so doing we get 

(-i^d u -im-^x)S = 5^(x-y). (20) 

Note that corrections to the interaction have the form of a series in powers of (l/(m + ax)T g ), and can be neglected 
in the first approximation. For this case both color-electric and color-magnetic terms contribute to the kernel M . 

Let us now turn to the second case, Eq. ([l9|). Since S(h' , x, y) varies as exp(— (ax + m)\h'\) (cf Eq.(p"l|)) and there 
is the factor exp(— (ax + m)\h — h'\) in (|16|), it is essential that the above factors being integrated out with 6(h — h') 
or 9(h' — h) in the expression ([l6]). For the term with 9(h — h') we get exp(— (ax + m)h) and bounds of integral are 
defined by T g , whereas for the term with 9(h' — h) the integration over dh' yields a factor l/(ax + m), so that this 



contribution does not grow for large x. As a consequence the color-magnetic term can be neglected in case (19) 
Consequently writing the S , (/i,x,y) in the form 

S(h, x, y) = ie-(-+ m ) W 5 (x, y) ( °W Q{ _ h) ) (21) 

where g(x, y) « <5' 3 ^(x — y) is a smeared <5-function, one obtains an equation 

{-ilud a -irn- iax)S(h,x.,y) = S ( - 4 \x - y), (22) 

where all interaction ax is due to the electric term J^ B ^(x,x). 

In this way we have confirmed a posteriori that solution S of (|l|) has at large distances the form of the Green's 
function for the linear potential, as given by (p2[), and hence our choice of S — Su n as the first approximation for the 
kernel M, Eq. (g), is justified. 

Let us now discuss the regimes ( |l8| ) and dl9] ) in more detail. Consider first the case of heavy quark mass, 

mT g > 1. (23) 

In this case one automatically obtains the regime ( JToj ) and hence the linear potential as in j22]). Since T g ~ lGeV -1 , 
only top and bottom quarks satisfy (|23), while charmed quark mass lies at the boundary. In the limit m — > oo the 
Green's function Su n Eqs. dqjll]) becomes the standard heavy-quark expression 

j p — m\h\ 

S lm -^S = S^(x - y) {6(h)(l + 74 ) + e(-h)(l - 74)} , (24) 

which agrees with Refs. |lj, @] and ||. However one should not interpret this as corresponding to an admixture 
of scalar and vector confining pieces, i.e. V m ix — ca;(| + 574), since as it was shown explicitly in |HJ, and also here, 
in view of the symmetry of the spectrum, that in the limiting case mT g 3> 1 the potential has to be pure scalar 

V sca lar = CTX, Eq. (||) . 

We now turn to the case (m + ax)T g « 1. It is clear that this condition is valid only for a restricted region of x. 
However, for light quarks the selfconsistent solution of (Q) with the replacement (0) in the kernel M is not a good 
approximation for calculating the spectrum of the lower lying states, since Eq. ( |23| ) is violated. This conclusion agrees 
with the one found in [l5| . For low mass states, where the effective region of interaction x e ff satisfies 

(m + ax eff )T g < 1, (25) 

one should use Eq. (|l^) with the approximation valid for the the regime mT g C 1. It is a reasonable starting 
approximation for the whole region of x as long as we are interested for the spectrum, satisfying the condition (|2q). 



It should be noted, that in that case the controversy discussed in Ref. jl5[ for the replacement (24) does not take 
place. Moreover, the conclusion reached in Ref. (lj] that the effective interaction has the form V = ^ax, applies only 
to the states for which x e ff satisfies Eq. (|25|), whereas for higher excited states inevitably another regime, Eq. ©, 
starts to apply with V e ff = ax instead of ^ax. 

The region of validity of the magnetic string tension is important from the physical point of view, since the additional 
\ax originates from magnetic field correlators. From our analysis we clearly see that at asymptotic large x, i.e. is 
for very long (and therefore heavy strings), the confining mechanism is purely electric. Moreover, the field contents 
is independent of the quark masses at the end of the string. It is only at intermediate regions where the magnetic 
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contribution may play an important role. In particular, we have found that instead of regimes mT g 3> l,mT g <C 1 
investigated in Refs. jll} and [ ^4|Jl5| one has the two regimes (p!s|) or ( |l9| ) where the total mass of the string plus quark 
mass enters, and the resulting confining force is linear, but with different strength. For heavy quarks with mT g 3> 1 
the regime (|l8| ) is essentially absent and we may safely use the color electric confinement mechanism. 

From the phenomenological point of view the lowest states of light quarks with the property ( |l5| ) feel string tension 
|<7 and this may be important for the resulting masses of heavy-light systems, such as D, D s , B, B s as it will be 
demonstrated in the next section (a similar remark about D s and B s was made earlier in [[L4|). 

Moreover, this increase of string tension resolves substantially (at least for the lowest levels) the descrepancy found 
for the Regge slopes of light mesons using relativistic quasi-potential equations for particles with spin, c.f. |17| ]. In 
particular, it was found in Ref. [jl6| that a larger string tension of a = 0.33 GeV 2 than the usually accepted value 
of 0.18 GeV 2 was needed to fit the experimental spectrum of the light mesons. This should be contrasted with the 
prediction ct\ — for the Regge slope, given by J = «o + «iM 2 , in the case of the spinless Salpeter equation p8| , 
which is close to the nonrelativistic prediction. A similar result is found in the case of our nonlocal kernel for the 
light-heavy quark system. Considering the Regge trajectories as obtained in Ref. |l3|] and Table I of this paper we 
find for the case of the light-heavy quark system a value of a\ = — . This is to be compared with the spinless Salpeter 
prediction for this system, given by a\ — 4-. The effect of the covariant treatment of the spin can already be seen 
when we use the linear Dirac equation with V(x) = ax + c for heavy-light mesons. In this case we get a slope of 
Y~ for Co = 0. Most of the difference between our results and those of the linear Dirac equation can be attributed 
to having at large distance an effective negative constant term cq in the linear potential present in the case of the 
nonlocal kernel |13| , which leads to the shifting of the bound state masses. It should be noted, that there is another 
mechanism ]l9| , |20[ to decrease the Regge slope. Considering a rotating string [|l9| the Regge slope gets somewhat 
closer to the nonrelativistic result. One finds a value of — in the case of the light-heavy quark system, corresponding 
to a physical half-string. 

III. SPECTRUM OF HEAVY-LIGHT MESONS 

The analysis of the previous section suggests that the lowest bound state solutions of (|l|) can be determined from 
the approximate instantaneous nonlocal Dirac equation of the form 

(ap + /3m)^W +P / M(x,z)?A n (z)d 3 z = £„V™(x). (26) 



Using the WKB solution for the Green's function S we find that for small T g the kernel M can be approximated by 



M(x,z) = 
where i5(x,z) is defined as 



V^T s j( B '(x,z)+ / dh'8(h')e ^e- 2 ^ +m ) ,l '^ fe M) (x,z) 7l7fe 



<5(x,z), (27) 



<5(x,z) = -^-*(l-co8d M )[Ao(o)+Ai(a)] (28) 

47TZ 



with Ao(a) and Ai(et) denoting the limiting values of Ao(a) and Ai(a) respectively when A = axh ~ axT g tends to 

zero. For low lying states of the light quark system the spatial regions of interest are expected to satisfy the condition 
(|l8|), i.e. (m + ax)T g <C 1. Hence we get 

M(x,z) = v^T 9 {J^)(x,z) + 4 M) (x,z) 747fe }5(x,z). (29) 

It is seen in ( p8| ) that at large x and z the function <5(x, z) tends to S^(x. — z) and from ( p7| ) one immediately 
obtains 

M(x, z) ^crxS {3) (x - z) (30) 



In what follows we shall solve Eq. (26) in two cases i) when M(x, z) is replaced by its large distance limit ( J30| ) ii) 
when M(x, z) is approximated by 
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M(x,z) = M l(X)Z ) = 0Fr4j g (x,z) ^ (1 2 ^^ Koia). (31) 




The latter approximation is justified in the situation when |x — z| <C x, which follows from the exponential damping 
of Kq(cl) when \a\ grows. 

One can exploit the computations from |H| to obtain the eigenvalues £ n {j,l) of Eq. (|2^) with the kernel for 
the lowest states, given in Table 1 for the case of the vanishing quark mass, m = 0, and in Table 2 for the mass 
to = 0.15 and in Table 3 for to = 0.20 GeV. 

In what follows we shall first concentrate on the s — and p— state eigenvalues and compare them to the results of 
the QCD sum rules and lattice calculations. In the language of heavy quark effective theory one has the following 
expansion ]22|- ]25| for the heavy-light meson mass ran 

m h = m Q \l + A- + -L-(X 1+dH X 2 ) + 0[ — ) ), (32) 

where A(n, j, /) = e n (j, I). Using a = 0.18 GeV 2 the solution of (|2|) with the kernel ( |3l| ) yields the S-wave eigenvalue 
A(0,±,0) = A s 

A s = 0.520 GeV for a s = 0, A s = 0.360 GeV for a s = 0.3, (33) 

which are about a factor of y| smaller in the absence of the magnetic contribution. The predicted values (|33| ) should 

be compared with the results of the QCD heavy-flavour sum rules [p2]-[25|l As = 0.57 ± 0.07 GeV and the result of 
the analysis from semileptonic B decays [^6| As = 0.39 ±0.11 GeV. A similar value was obtained recently from the 
QCD sum rules @- 

As = 0.45 ±0.15 GeV (34) 

One can see a reasonable agreement of our results with the latest sum rule calculations ( |34| ) . Note here that we have 
taken into account the color Coulomb interaction to all orders, whereas in the QCD sum rules only the leading order 
term is retained, therefore one may expect that higher orders will decrease somewhat the value (|3~4|). 

In a similar way one may compute energy eigenvalues for the strange heavy-light mesons. From Tables 2 and 3 we 
see that with a = 0.3 we have 

A^ s) = 0.445 GeV for to = 0.15 GeV, A£° = 0.476 GeV for to = 0.20 GeV. (35) 

These numbers can be compared to the values from the experimental B s and D s masses. We find 

AM S (B) = M Bs - M B = 90 MeV, AM S (I5) = M Ds - M D = 100 McV. (36) 

Similar values are found from the spectrum of heavy-light mesons, computed recently on the lattice [ p8| for strange 

mesons. From Eq. ([55|) we see that the experimental data are close to our predicted value of AM S — Kg — As = 
85 MeV for to = 0.15 GeV. The various available data on As are summarized in Table 4. 

We turn now to orbital and radial excitations. For the states with 1 = 1, and j = | and \ the mass splitting is due 
to the spin-orbit interaction inherent in the Dirac equation. Denoting these energies as s n (j, 1) = Ap(j) we find for 
the nonstrange quark (in GeV): 

A P (i) = 0.817, A P (|) = 0.732 {a, = 0) (37) 

Ap(~) = 0.665, A P (|) = 0.620 (a a = 0.3). (38) 

Similarly for strange mesons with a strange quark mass to = 0.15 GeV, we obtain 

A^ s) (i) = 0.898, A^(|) = 0.832 (a s = 0) (39) 

A^(~) = 0.741, Ap(|) = 0.712 (a s = 0.3). (40) 

These calculations can be compared with the results of lattice calculations in |28| , with experiment and the recent 
QCD sum rule calculations p7|.The latter yield for m = 
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Ap = (1 ± 0.2) GeV. 



(41) 



This value is somewhat higher than the results (37) and (pS)). Lattice calculations in [28] give for the difference 
M(Bj) - M(B) ~ Ap - A s « 456 MeV, which should be compared with our results, AA = A P (±) - A s ~ 305 MeV 
for a s = 0.3 and with experiment M(B}) ~ M{B) ~ 338 MeV. Here M(B) = \M B (\-) + \M b (Q~). In addition 



there is a calculation of heavy-light mesons in the framework of the QCD string approach |29|, where the only input is 
current quark masses (m„, m,i, m s ), string tension a and a s . These results have been obtained in ]30| ] and recently in 
ifjlf for real B, B s , D, D s mesons and are easily computed for the limiting case of m q — > oo, which yields values listed 
in Table 5. The rather low value found for A A suggests that we still miss some strength in the orbital excitation in 
the present work. In Table 5 the results of different approaches to Ap and Ap- 1 are collected. 

Radial excitations are readily obtained from solving Eqs. ((2^), (|3l]) and yield for the n = 1 state 

£i(~,0) = 0.951 (a, = 0), 0.805 (a, = 0.3) (42) 
and for the strange meson with m = 150 MeV 

= L036 ( a * = 0), 0.880 (a s = 0.3) (43) 
while for the radial excitation with I = 1 one obtains 

£l (j, 1) = 1.140 (a s = 0), 0.997 (a s = 0.3) (44) 
e[ s) (j, 1) = 1.221 (a, - 0), 1.076 (a, = 0.3) (45) 

These values are compared in Table 6 with the results of the QCD string approach and recent lattice calculations 



IV. SUMMARY AND CONCLUSIONS 



The main results of the present paper are of both theoretical and phenomenological scope. On the theoretical side 
it is shown in Section 2, that there exist two possible dynamical regimes for the quark at the end of the Dirac string, 
namely (18) and (19). For the case of heavy quark with rnT g <C 1 only one regime (19) is available and results are 
the same as discussed in [11,14,15]; i.e. color magnetic contribution to the string is suppressed in this case and one 
has at large x local Dirac equation with linear potential. 

For the case of light quark there is a possibility of another regime (18), where color magnetic field also contributes. 
It is shown that this regime operates at intermediate distances and yields a static Dirac equation with increased string 
tension. It is demonstrated also that the use in this case in the kernel of nonlinear Dirac equation of the full quark 
propagator S (or its WKB approximation Swkb) leads to the consistent results, while the use of the heavy-mass 
propagator So - leads to inconsistencies shown in [15]. On the phenomenological side the regime (18) yields the energy 
eigenvalues which are in a better agreement with other calculations and experiment, as demonstrated in Tables 4-6, 
as compared to our previous results [13] where color-magnetic contribution has been neglected. 

The study of the role of color magnetic fields in the dynamics of light quarks is at its beginning and the first results 
call for more detailed investigation of the transition between regimes (18) and (19) and other applications, e.g. to 
mesons and baryons consisting of light and heavy quarks. 

One of the authors (Yu.S) is acknowledging a partial financial support through the RFFI grants 00-02-17836 and 
00-15-96785. He is also grateful to Yu.S.Kalashnikova and A.V.Nefediev for useful discussions and the use of their 
data prior to publication. 
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Table 1 





io 




1-1 


2 2 

2> z 


2 2 

2' z 


2 3 


i 3 

9 ! d 


a s = 
a s = 0.3 


0.520 
0.360 


0.817 
0.665 


0.732 
0.620 


0.934 
0.885 


0.911 
0.818 


1.147 
1.057 


1.070 
0.987 



Ground state energy eigenvalues s n (j, I) in GeV, n = 0, for two values of a s and cr = 0.18 GeV 2 , T g = 0.23 fm, m=0. 

Table 2 





£,0 




2 1 


2 2 

2> z 


2 2 

2' z 


5 Q 
9 j ° 


i 3 


a s = 
a 8 = 0.3 


0.623 
0.445 


0.898 
0.741 


0.832 
0.712 


1.078 
0.965 


1.010 
0.911 


1.233 
1.140 


1.168 
1.081 



The same as in Table 1 but for m = 0.15 GeV. 

Table 3 





i,0 




2 1 

2' 1 


2 2 

2> Z 


2 2 

2' Z 


2 3 

2'° 


1 3 

2 ' J 


a s = 
a s = 0.3 


0.659 
0.476 


0.927 
0.769 


0.867 
0.744 


1.107 
0.994 


1.044 
0.944 


1.263 
1.169 


1.202 
1.114 



The same as in Table 1 but for m = 0.20 GeV. 
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Table 4 



Energy eigenvalues A5 of the heavy-light system in the static heavy quark approximation obtained in different 

approaches. 



Refs. 


Method 


A s (GeV) 


24 


QCD sum rules 


0.5 


25 


QCD sum rules 


0.4-7- 0.5 


27 


QCD sum rules 


0.45± 0.15 


26 


Experiment 


0.39± 0.11 


13 


Nonlin. Dirac 


0.287 


this work 


Nonlin.+ magnetic 


0.360 



Table 5 



The same as in Table 4 but for Ap — 



Refs. 


Method 


A P - A 5 (GeV) 


27 


QCD sum rules 


0.55 ± 0.35 


31 


QCD string 


0.40 


28 


Lattice 


0.47 


PDG 


Experiment 


0.383 


this work 


Nonlin. + magnetic 


0.305 (Ap(i)-A s ) 



Table 6 

The same as in Table 4, but for the radial excitation, A' s — Ag. 



Refs. 


Method 


A' s - A s (GeV) 


32 


experiment, Mb* — Mb 


0.581 


28 


lattice, Mb* — M B 


0.602 


31 


QCD string 


0.564 


this work 


Nonlin. +magn. 


0.631 
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